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C3 . ABSTRACT 

^ '. 

. £^ | Aims. We present the first large, unbiased sample of Lyman Break Galaxies (LBGs) at z ~ 1. Far ultraviolet-dropout (1530 A) 
i galaxies in the Chandra Deep Field South have been selected using GALEX data. This first large sample in the z ~ 1 universe 
5_l ' provides us with a high quality reference sample of LBGs. 

\ Methods. We analyzed the sample from the UV to the IR using GALEX, SPITZER, ESO and HST data. 
Results. The morphology (obtained from GOODS data) of 75 % of our LBGs is consistent with a disk. The vast majority of 
LBGs with an IR detection are also Luminous Infrared Galaxies (LIRGs). As a class, the galaxies not detected at 24pm are 
an order of magnitude fainter relative to the UV compared with those detected individually, suggesting that there may be two 
types of behavior within the sample. For the IR-bright galaxies, there is an apparent upper limit for the UV dust attenuation 
and this upper limit is anti-correlated with the observed UV luminosity. Previous estimates of dust attenuations based on the 
ultraviolet slope are compared to new ones based on the FIR/UV ratio (for LBGs detected at 24 /jm), which is usually a more 
reliable estimator. Depending on the calibration we use to estimate the total IR luminosity, /3-based attenuations Afuv are 
larger by 0.2 to 0.6 mag. than the ones estimated from FIR/UV ratio. Finally, for IR-bright LBGs, median estimated /3-based 
SFRs are 2-3 times larger than the total SFRs estimated as SFRtot = SFRuv + SFRir while IR-based SFRs provide values 
below SFRtot by 15 - 20 %. We use a stacking method to statistically constrain the 24/xm flux of LBGs non individually 
detected. The results suggest that these LBGs do not contain large amounts of dust. 
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1. Introduction 
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Lyman Break Galaxies (LBGs) are the most numerous ob- 
jects observed at high redshift (z > 2 — 3) in the rest-frame 
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ultraviolet (UV). The discovery of a large population of 
LBGs beginning with the work of Madau et al. (1996), 
followed by the spectral confirmation of their redshifts, 
provided the astronomical community with the first large 
sample of confirmed high redshift galaxies (Steidel et al. 
1996; Lowenthal et al. 1997). The spectra of bright LBGs 
(e.g. cB58 by Pettini et al. 2000; Teplitz et al. 2000) are 
remarkably similar to those of local starbursts, indicat- 
ing that these objects are forming stars at a high rate. 
The observed colors of LBGs are redder than expected 
for dust-free star-forming objects. This reddening suggests 
that some dust is present in this population. However, the 
amount of dust in LBGs (Baker et al. 2001; Chapman 
et al. 2000), and therefore the reddening-corrrected star 
formation rate (SFR C ), is poorly known. Meurer et al. 
(1999), Adelberger & Steidel (2000) and subsequent pa- 
pers tried to estimate the amount of dust attenuation from 
the j3 method (Calzetti et al. 1994). However, it has been 
shown recently that this approach only provides rough es- 
timates of the total UV attenuation in local galaxies (e.g. 
Buat et al. 2005; Burgarella et al. 2005; Seibert et al. 2005, 
Goldader et al. 2002 and Bell 2002). 

High redshift LBGs are mainly undetected at sub- 
millimeter (sub-mm) wavelengths, where the emission of 
galaxies is dominated by the dust heated by young stars 
(Kennicutt 1998). Only the most extinguished LBGs are 
detected by SCUBA (Chapman et al. 2000, Ivison et al. 
2005) and we have no idea of the dust attenuation for a 
representative sample. Very recently, Huang et al. (2005) 
observed a population of LBGs at 2 < z < 3 detected 
with SPITZER. Unfortunately, due to the very high red- 
shift, the SPITZER/MIPS observations were not deep 
enough to detect the thermally reradiated emission from 
very many LBGs at z ~ 3 and the SPITZER/IRAC 
data, although deep enough to detect many LBGs, only 
probe the rest-frame NIR (i.e. no information about the 
dust enshrouded star formation can be inferred) . The mor- 
phology of LBGs is also a matter of debate: early works 
(e.g. Giavalisco et al. 1996) suggested that LBGs could 
be ellipsoidal, and therefore perhaps the progenitors of el- 
lipticals or of the bulges of spiral galaxies. The problem 
is that we can hardly detect low surface brightness areas 
at high redshift because of the cosmological dimming. For 
instance, Burgarella et al. (2001) suggest that only com- 
pact star forming regions could be easily detected in deep 
HST observations. 

On the other hand, observations in the sub-mm range 
have revealed a population of FIR-bright galaxies that 
might be similar to local (Ultra) Luminous IR galaxies 
((U)LIRGs) (Blain et al. 1999) with 1O 11 L < L IR = 
L(8-1000/zm) < 10 12 i Q for LIRGs and 10 12 L Q < L IR < 
10 13 Lq for ULIRGs. These objects are likely to dominate 
the Cosmic Infrared Background (CIB; Elbaz & Cesarsky 
2003) at high redshift. The link between LBGs and LIRGs 
is still an open question: are there two classes of objects 
or are they related? If they are two facets of the same 
population, then we could, for instance, correct UV fluxes 



for the dust attenuation to recover the full Star Formation 
Density (e.g. Perez-Gonzalez et al. 2005). 

The usual way of detecting and identifying LBGs is 
through the so-called dropout technique, i.e. the absence 
of emission in the bluest of a series of bands due to the 
Lyman break feature moving redwards with the redshift 
(e.g. Steidel & Hamilton 1993, Giavalisco 2002). However, 
it is well known that selection effects can have a very 
strong influence on the deduced characterics of an ob- 
served galaxy sample (e.g. Buat et al. 2005; Burgarella 
et al. 2005). Until now, there has been no way to detect 
a general, unbiased sample of LBGs at low redshift (i.e. 
z < 2), with the same dropout method very successfully 
used at z > 2 because we lacked an efficient observing 
facility in the UV range. This was quite unfortunate be- 
cause the simple fact that the galaxies are closer to us 
means that we can access much more information on the 
morphology, detect fainter LBGs in the UV and in the 
IR, and therefore harvest larger samples. GALEX and 
SPITZER changed this situation and have allowed us to 
define a large sample of LBGs at z ~ 1 in the Chandra 
Deep Field South (CDFS). 

In this study, we combine the detection in the UV of 
true (i.e. with a detected Lyman break) LBGs at z ~ 1 
with GALEX and at 24/zm with SPITZER/MIPS. These 
data let us estimate the total dust emission and therefore, 
the dust-to-UV flux ratio, which provides a good tracer 
of the dust attenuation in the UV. We also use high spa- 
tial resolution images to analyse their morphology. We are 
therefore able to perform a complete analysis for the first 
time on a large LBG sample. 

A cosmology with H = 70 km.s -1 .Mpc -1 , = 0.3 
and Qa = 0.7 is assumed in this paper. 

2. The Galaxy Sample 

GALEX (Martin et al. 2005) observed the CDFS field 
for 44668 sec (Deep Imaging Survey = DIS) in both the 
far ultraviolet (FUV) and the near ultraviolet (NUV). 
The GALEX field is centered at a = 03h32m30.7s, <5 
= -27deg52T6.9" (J2000.0). The GALEX IR1.1 pipeline 
identified 34610 objects within the 1.25°-diameter field of 
view. The GALEX resolution (image full width at half 
Maximum = FWHM) is about 4.5 arcsec in FUV and 6 
arcsec in NUV. 

This field has also been observed by SPITZER using 
MIPS (Rieke et al. 2004) in the guaranteed time observ- 
ing program. The MIPS observations provide about 7-8 
sources arcmin~ 2 at 24 [im centered on a ~ 1.45 x 0.4 = 
0.6 deg 2 field of view. The SPITZER image FWHM is 
about 6 arcsec and almost perfectly matches that of 
GALEX. 

Redshifts from GOODS (Vanzella et al. 2005) and 
VVDS (Le Fevre et al. 2004) are available at the cen- 
ter of the GALEX field. Part of the GALEX CDFS 
field was observed by COMBO 17 (Wolf et al. 2004) 
over 0.5 x 0.5 deg 2 . We made use of COMBO 17 red- 
shifts for objects with r < 24.5. In this range, the quality 
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Sz/(1 + z) remains within a z < 0.03 for 53 % of the ob- 
jects. Finally, we obtained photometry from the European 
Southern Observatory Imaging Survey (EIS) in U, B, V, 
R and I. 

We built a sample of LBGs as follows. From the sources 
with redshifts, we selected objects in the rest-frame FUV 
(i.e. in the observed NUV) that we cross-correlated (r = 1 
arcsec) first with the ground-based optical EIS data, then 
(r = 4 arcsecs) with the MIPS 24 fim data. In the resulting 
catalogue, we down-selected to objects that COMBO 17 
puts in the "GALAXY" class. Then, we extracted the 
sources with redshift 0.9 < z < 1.3. Since we wish to 
study Lyman break galaxies, we omitted objects without 
observed NUV and U flux : we kept only objects down 
to GALEX NUV magnitudes = 24.5 corresponding to the 
GALEX 80% completeness level and the U-band limiting 
magnitude at U=25.1. We did not use a color-color se- 
lection as performed by Steidel & Hamilton (1993). Since 
we are studying galaxies with known redshifts, the extra 
color is not needed to screen out interlopers. The selec- 
tion on the x-axis (i.e. G — R color) might bias the sample 
toward low-reddening LBGs that we wish to avoid so we 
can achieve a more general understanding of the sample. 
In fact, we find that the members of our sample fall within 
the traditional color-color LBG range, or very close to it, 
as discussed in Section 3.5. For galaxies at z ~ 1, GALEX 
FUV corresponds to a rest-frame wavelength of ~ 765 
A and GALEX NUV corresponds to - 1155 A. The ob- 
served FUV and NUV filters are therefore in the same rest- 
frame wavelength ranges as U and G filters used to identify 
Lyman breaks at high redshift (e.g. Giavalisco 2002 for a 
review). The observed FUV-NUV color thus gives a clear 
indication of the Lyman break. We picked the objects with 
the strongest indication of a break : FUV ~ NUV > 2; 
the final sample contains 297 LBGs. Of this list, 49 ob- 
jects (16.5 %) have a measured flux above the 80 % com- 
pleteness limit of the SPITZER 24 \xm data (e.g. Perez- 
Gonzalez et al. 2005) and thus at a high enough ratio of 
signal to noise to be treated individually. Fig. 1 shows two 
examples of these LBGs at several wavelengths including 
the two GALEX bands, the B and R EIS bands and the 
24/xm SPITZER/MIPS band. Many additional galax- 
ies were detected by SPITZER but at a weaker level; 
below we will describe how we used a stacking technique 
to probe their properties. Our sample of 297 UV-selected 
LBGs in the redshift range 0.9 < z < 1.3 constitutes 
the database that we will study in this work (except for 
Sect. 3.5). The common GALEX - SPITZER/MIPS - 
COMBO 17 field of view (mainly limited by COMBO 17) 
is about 0.25 deg 2 , which translates to ~ 1180 LBGs deg -2 
and ~ 200 LBGs deg -2 for which we have individual IR 
detections. 

3. Lyman Break Galaxies at z ~ 1 

Our LBGs provide the first opportunity to study an un- 
biased sample of LBGs systematically at z ~ 1. We ana- 
lyze the UV and IR luminosities of these galaxies, measure 
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Fig. 2. The distribution in observed luminosity presented 
here corresponds to UV luminosities in blue (solid) and IR 
luminosities in red (dashed). Heavy lines are drawn from 
the population of LBGs with detected counterparts in the 
24/zm MIPS image and thin lines to upper limits only. 
Note that the cut at low Lju is not sharp because the 
83/iJy limit used here is not the detection limit but the 
80 % completeness limit. UV luminosities cover the same 
range for the two samples and reach uncorrected UV lumi- 
nosities of LogLjjy — ll(mLQ). This upper limit is con- 
sistent with Adelberger & Steidel range for LBGs at z ~ 3 
but seems inconsistent with the Balmer-break sample at 
z ~ 1. Two areas are shaded. Starting from the left, the 
first corresponds to the range covered by UV Luminous 
Galaxies (open right ended) and the second one corre- 
sponds to Luminous IR Galaxies. A few Ultra Luminous 
IR Galaxies are also detected. 



their morphologies and their star formation as revealed by 
the UV and IR data and finally discuss the implications 
they have for studies centered on higher redshifts. 

3.1. Ultraviolet and Infrared Luminosities 

We find a large range of observed (i.e. un-corrected for 
dust attenuation) FUV luminosities (A/a in rest-frame 
FUV) with 9.3 < LogL uv [L Q ] < 11.0. The lowest lu- 
minosity is set by the limiting magnitude in the U-band 
at U = 25.1. Below the break, the limiting magnitude 
amounts to FUV — 26.0. Although fainter objects are 
detected in the NUV (down to NUV = 25.9), we use 
a limiting magnitude of NUV = 24.5 to compute safer 
FUV — NUV colors and therefore perform a safer Lyman 
Break selection. The average value is < LogLjjv[LQ] > = 
10.2±0.3 for the sample with individual IR detections, and 
< LogLjjv[L©] >= 10-1 ± 0.3 for the rest of the sample. 

Total IR luminosities {Lib) are estimated following 
the procedure described in Perez-Gonzalez et al. (2005). 
Briefly, rest-frame 12 /im fluxes are calculated by com- 
paring the observed mid-IR SEDs (including IRAC and 
MIPS fluxes) of each individual galaxy with models of 
dust emission (e.g. Chary & Elbaz 2001 or Dale et al. 
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Fig. 1. Two galaxies from our LBG sample are shown here, from left to right in GAL EX FUV and NUV, then EIS 
B, HST GOODS B, EIS I, HST GOODS I and finally SPITZER/MIPS 24^m band. For the two galaxies, the 
leftmost image (blue frame) is below the Lyman break at z ~ 1 and the galaxy is not visible. The size of GOODS 
images is 4 x 4 arcsec 2 . The corresponding GOODS field (green frame) is plotted in the large (35 x 35 arcsec 2 ) 
EIS images, a) a LBG classified as a disk-dominated galaxy; the more compact object is at z = 0.546 from VVDS, 
it appears reddish and should not contribute in ultraviolet and in far infrared; b) a LBG classified as a merger / 
interacting galaxy. 



2002). This procedure is meant to cope with the strong 
K-corrections observed in the mid-IR due to the emis- 
sion from aromatic molecules. We use the formulations of 
Takeuchi et al. (2005) and Chary & Elbaz (2001) for con- 
version from 12y«m flux density to Lir. In addition to the 
intrinsic differences due to the two calibrations (the for- 
mer provides Lir lower by 0.2 dex), the conversion from 
L\2^,m to Lm can introduce errors up to a factor of 2 for 
individual normal galaxies and 4 for galaxies with SED 
variations over the full IRAS sample range (Takeuchi et 
al. 2005; Dale et al. 2005). The effects of these errors are 
greatly reduced in this work because we discuss average 
properties, not those of individual galaxies; in this case, 
the uncertainties are likely to be commensurate with those 
in the overall conversions to Lir (i.e., still ~ 0.2 dex). 
There are also uncertainties in the luminosity values due 
to the distance (since we use photometric redshifts with 
6z/(z + 1) < 0.03), but they are less than 10% for Ljjy 
and Lir. 

It is difficult to compare our sample to previously pub- 
lished ones since none was available in this redshift range 
before GALEX. However, Adelberger & Steidel (2000) 
have discussed a Balmer-break sample at z ~ 1. The dis- 
tributions in Ljjv and Lir for our sample are shown in 
Fig. 2. The lower limit is set by the flux limits but the 
upper limit of our LBG sample is about the same as the 
z > 3 one in Adelberger & Steidel's (2000). However, the 
upper limit of our sample is higher by a factor of ~ 4 
(assuming H = 70) than Adelberger & Steidel's (2000) 
sample. 

Heckman et al. (2005) defined UV Luminous Galaxies 
(UVLGs) as galaxies with UV luminosities above 
Log(Luv) = 10.3 Lq. They found that these UVLGs bear 
similarities to LBGs, especially a sub-sample of compact 
ones. In our sample, 22.2 % of the LBGs are UVLGs, 30.6 
% of the LBGs with an IR counterpart are UVLGs. The 
83/iJy detection limit at 24/mi approximately corresponds 



to Log(Lm) > IILq, that is, nearly all the IR-detected 
LBGs are LIRGs (95.9 %) and there are 2 ULIRGs (4.1 
%). An association between UVLGs and LIRGs was pro- 
posed by Burgarella et al. (2005). We confirm here this 
association and extend it to LBGs. 

3.2. The Morphology 

The Great Observatories Origins Deep Survey (GOODS) 
provides high resolution and high signal-to-noise images 
of some of our LBGs, which can be used to study their 
morphology in the rest-frame B band. An advantage of 
our low redshift sample of LBGs is that the images ex- 
tend to low surface brightness and hence morphologies 
can be determined well. We compute the asymmetry and 
concentration (Fig. 3) as in Lauger et al. (2005a) from 
the objects within the GOODS field which have a signal- 
to-noise ratio larger than S/N « 1 per pixel and whose 
coordinates are within 2 arcsecs from the GALEX detec- 
tion. We were able to obtain the morphology for only 36 
LBGS out of our 300 LBGs (about 1/4 of our GALEX 
+ SPITZER + COMBO 17 field is covered by GOODS). 
Fig. 3 leads us to two conclusions: i) all but one LBG in 
our sample are located on the disk side of the line sepa- 
rating disk-dominated and bulge-dominated galaxies, and 
ii) part of them (22 %) are in the top part of the dia- 
gram, i.e. with an Asymmetry larger than 0.25 and could 
be interpreted as mergers. This kind of quantitative anal- 
ysis is also applied to higher redshift LBGs, however, we 
must be careful in the interpretation because, even if disks 
are present, it would be very difficult to detect them due 
to the cosmological dimming (e.g. Burgarella et al. 2001). 
Indeed, in deep spectroscopic observations (e.g. Moorwood 
et al. 2000, Pettini et al. 2001), the profiles of the optical 
nebular lines suggest the presence of disks in some LBGs. 

A number of studies have been devoted to galaxy mor- 
phology in the redshift range 0.6 < z < 1.2. At z ~ 0.7, 
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most of the works seem to agree that about 60 - 70% 
of the objects can be classified as disk-dominated galaxies 
(mainly spirals and Magellanic irregulars) and 10 - 20% as 
mergers/interacting galaxies. Here spirals are a sub-group 
of disks which exhibit a more symmetric (spiral) struc- 
ture than irregular-like objects similar to the Magellanic 
clouds. Their asymmetry is therefore lower. Lauger et al. 
(2005b) found about 70 - 80% of disk-dominated galaxies 
at z ~ 1. Zheng et al. (2004) studied the HST morphology 
of a sample of LIRGs and also found that a large major- 
ity (~ 85 %) of them are associated with disk- dominated 
galaxies. This conclusion is reached whether the selection 
is in the rest-frame ultraviolet (Wolf et al. 2005) or in the 
infrared (Bell et al. 2005). However, some dispersion due 
to the cosmic variance might exist (Conselice, Blackburne 
& Papovich 2005). 

Therefore, overall it appears that the majority of the 
star formation at 0.6 < z < 1 resides in disks and about 
half of it in spirals. The numbers that we draw for our 
LBG sample at z ~ 1 are globally consistent: we find 
that ~ 22% of the LBGs are likely mergers (e.g. Fig. lb), 
~ 75% are disks (e.g. Fig. la) and only ~ 3% (i.e. 1 
galaxy) is possibly a spheroid. In the cases where our 
LBGs can also be classified as LIRGS from their IR lumi- 
nosity, our measured morphologies are consistent with the 
LIRG morphology measurements of Zheng et al. (2004). 

3.3. Ultraviolet Dust Attenuations 

Until now, it has been difficult to estimate the validity 
of dust attenuation estimates for distant LBGs, because 
we had no clear idea of their Lir. Adelberger & Steidel 
(2000) tried to estimate the 800/im fluxes of their LBG 
sample from the (3 method and compared the results to 
observations. However, only the most extreme LBGs can 
be detected either directly in the sub-millimeter range or 
in the radio range at 1.4 GHz and therefore could be used 
in this comparison. 

In this paper, we use total IR luminosities, £/_r, and 
Luv to compute the FIR/UV ratio, which is calibrated 
into FUV dust attenuation Apjjv (e.g. Burgarella et al. 
2005). This method has been shown to provide more ac- 
curate dust attenuations than those from the UV slope (3. 
Fig. 4 shows an apparent anti-correlation of Apuv with 
the UV luminosity. It is not clear, however, whether this 
relationship is real or only observational. Indeed, in ad- 
dition to the observational cut at low LpjjVi the 24//m 
lower limiting flux means that we cannot detect individ- 
ually low-luminosity galaxies with low dust attenuations. 
It is very interesting to note that we do not detect LBGs 
with both a high UV luminosity and a high UV dust atten- 
uation, and this cannot be caused by observational limits. 
In other words, we seem to observe a population of high 
Luv LBGs (which qualify as UVLGs) with dust attenua- 
tions similar to UV-selected galaxies in the local universe 
(e.g. Buat et al. 2005). UVLG galaxies are LIRGs with 
the lowest Apuv- 
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Fig. 3. The morphology of the LBG sample is quantita- 
tively estimated from the asymmetry and the concentra- 
tion (sub-sample drawn from a larger CDFS analysis by 
Lauger et al. 2005b). The smaller number of objects pre- 
sented here is due to the smaller field of view of GOODS 
as compared to ours. Black circless are LBGs without 
IR counterparts while red diamonds represent LBGs with 
a 24/im MIPS detection. The line corresponding to the 
limit between disk-dominated and bulge-dominated galax- 
ies (Lauger et al. 2005a). The location in the diagram re- 
flects the morphological type of the galaxies: more asym- 
metrical LBGs (e.g. mergers) are in the top part of the 
diagram (A > 0.25) while early-type spirals would have 
A < 0.1. The LBG sample is mainly dominated (75 %) by 
disk-dominated galaxies and the contribution from merg- 
ers amounts to ~ 21 %. 



There are now studies of high-redshift LBGs (z > 2) 
with SPITZER (e.g. Labbe et al. 2005; Huang et al. 
2005). However, the results are so far inconsistent: Labbe 
et al. (2005) found that LBGs are consistent with low- 
reddening models while Huang et al. (2005) found more 
reddened LBGs. Further observations will resolve these 
differences and provide a firm basis for comparison with 
our sample. 

With the 24^to SPITZER flux for the individually 
detected z ~ 1 objects, we can go a step further and check 
how UV dust attenuation estimations carried out from 
the (3 method compare with the better IR/UV-hased es- 
timates. This comparison has already been performed for 
nearby galaxies (Buat et al. 2005, Burgarella et al. 2005, 
Seibert et al. 2005 and references therein). Given the wide 
use of the (3 method on high redshift LBGs, it is useful to 
compare with our lower-redshift sample. 

Using the equations in Adelberger & Steidel (2000) 
(deduced from Meurer et al. (1999)), we estimate the IR 
luminosity that is used to compute Apuv and the total 
luminosity for each LBG. We observe a small overesti- 
mation of the dust attenuation as compared to the ones 
estimated from SPITZER/MIPS data and the dust-to- 
UV flux ratio. The /3-based mean dust attenuation esti- 
mated for our z ~ 1 LBG sample is Apuv — 2.76 ± 0.13 
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Fig. 4. Blue and red symbols are the same objects but 
luminosities are Luv for the former and Lj R for the lat- 
ter. The dust attenuation strongly decreases while Luv 
increases. Part of this apparent correlation might be due 
to the fact that observational limits prevent us from de- 
tecting low-luminosity LBGs with low dust attenuation. 
The clear cut on the upper parts of the box cloud cannot 
be due to observational biases. We do not seem to ob- 
serve UVLGs with high dust attenuations. On the other 
hand, the more dispersed but well-known increase of the 
dust attenuation with Li R is observed here. Larger sym- 
bols correspond to UVLGs (LogLjjv > 10.3Lq). Most of 
them have low dust attenuations but one is a ULIRG and 
has Afuv ~ 4. 



(er = 1.02) while the dust-to-UV dust attenuation gives 
A FU v = 2.16 ± 0.11 (a = 0.84) with L IR from Takeuchi 
et al. (2005) and A FUV = 2.53±0.12 {a = 0.94) with L IR 
from Chary & Elbaz (2001). The average value of the two 
dust-to-UV estimates is consistent with Takeuchi, Buat & 
Burgarella (2005). The net effect is that total luminosities 
based on the f3 method are slightly larger than the actual 
values and the deduced SFRs are therefore overestimated 
(see next section). The mean of the ratios of the (3 vs. 
FIR/UV A FUV = 1.31 ± 0.07 (a = 0.52) for Takeuchi et 
al. (2005) and A FUV = 1.09 ± 0.05 [a = 0.40) for Chary 
& Elbaz's calibration. 

Applying the Kaplan-Meier estimator (and using 
Chary & Elbaz's calibrations), we can take upper lim- 
its into account to estimate mean dust attenuations. 
We find moderate values < Apuv >— 1-36 ± 0.07 for 
Mpuv < —22 (5 data points of which 1 is an upper 
limit) and < A FUV >= 1.08 ± 0.11 for Mpuv < -21 
i.e. at z=3 (35 data points of which 65 % are upper 
limits). Using only detections, we reached, respectivelly, 
< Afuv >= 1-54 ± 0.09 and < A FUV >= 1.59 ± 0.18. 
Since we only use a small number of bright LBGs, this is 
hardly comparable to the numbers quoted in the previous 
paragraph. However, it suggests that lower Lpuv LBGs 
have higher dust attenuations (see Fig. 4). 



3.4. Star Formation Rates and Implications for the 
Cosmic Star Formation Density 

We estimate SFRs for our LBG sample from the IR lu- 
minosities and after applying dust corrections estimated 
from (3 and we compare them, in Fig. 5, to the total SFR: 
SFRtot = SFR UV + SFR IR where SFR UV is not cor- 
rected for dust attenuation. SFRtot is assumed to be 
the best SFR estimate and we use it as a reference. The 
first conclusion is that the dispersion is much larger for 
UV SFR c s computed with f3 dust corrections than for IR 
SFRs. But the median values are also different: SFR UV = 
112.3 ± 33.8 [a = 260.9) MQ.yr' 1 while SFR IR = 49.9 
± 13.1 (a = 100.9) MQ.yr' 1 if we use Chary & Elbaz 
(2001) and SFR IR = 30.5 ± 6.0 (a = 46.6) MQ.yr^ 
if we use Takeuchi, Buat & Burgarella (2005). Median 
SFR F o F for the two above calibrations are, respectiv- 
elly, SFRtot = 62.8 ± 13.2 (a = 102.1) MQ.yr' 1 and 
SFRtot = 41.1 ± 6.3 (a = 48.6) MQ.yr' 1 . As expected, 
for LBGs detected in IR, SFRj R is therefore a better es- 
timate. About 22 % of our LBGs with an IR detection 
have SFRtot > 100 MQ.yr' 1 (using Chary & Elbaz's 
calibration) as compared to less than 1 % in Flores et al.'s 
(1999) galaxy sample which confirms that our LBGs are 
forming stars very actively. However, none of the LBGs 
undetected at 24^im is above SFRtot > 100 MQ.yr' 1 . 

The higher SFRs reached when dust attenuations are 
computed with the UV slope j3 (depending on the Lj R 
calibration, +79 to +173%) lead to an overestimated con- 
tribution of LBGs to the Cosmic Star Formation Density if 
the same quantitative difference exists at higher redshift. 
However, Takeuchi, Buat & Burgarella (2005) showed that 
the current assumption of a constant dust attenuation 
docs not seem to be verified. The increase of the mean 
Afuv from 1.3 to 2.3 from z — to z = 1 means 
that, for a given observed FUV luminosity density, the 
dust-corrected star formation density would vary. Note 
that those mean dust attenuations cannot be compared 
with the numbers given for the Kaplan-Meier estimates 
which are biased toward large Mpuv LBGs while fainter 
LBGs seem to have larger dust attenuations in our sam- 
ple. Although FIR data are not always available, it is very 
important that one is aware of these uncertainties when 
using Star Formation Densities derived from UV values 
corrected from the f3 method for LBGs with high dust 
attenuations, especially at high redshift where we have a 
very poor knowledge of actual attenuations. These ambi- 
guities may be reduced by further study of Spitzer data, 
and with Herschel. 

3.5. Extension to Galaxies Faint at 24fim 

So far, most of our arguments have been based on the 
49 galaxies individually detected at 24/mi well above the 
completeness limit. To put the behavior of these galaxies 
in a broader context, we have determined average infrared 
flux densities for groups of galaxies by stacking. Images at 
24/im are shifted to a common center on the basis of the 
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SFR C „ V & SFRj 



Fig. 5. Blue open circles compare SFR^ V to SFRtot 
while red filled boxes compare SFRfir to SFRtot- Both 
SFRs are computed from Kennicutt (1998). The median 
SFRfir is underestimated by about 80 % and the un- 
derestimation increases at lower SFRs, which is consistent 
with the fact that the UV contribution (not accounted for 
from Lir) is usually higher at low SFRs than at high 
SFRs. For this sample of LBG detected at 24 /xm, the 
median value of SFRfj V is overestimated by a factor of 
2 - 3 in average with a possible trend for the difference 
to increase at high SFRs. The dispersion of /3-based UV 
SFR c s {usfr ~ 200M Q .yr- 1 ) is much larger than IR- 
based SFR (ctsfr ~ 30 — 40MQ.yr~ 1 depending on the 
calibration into Lm). 

24/im coordinates if the object was well enough detected 
in that band, or the coordinates of the optical identifi- 
cation otherwise. We use sigma-clipping to eliminate sur- 
rounding sources; the level at which clipping occurs is ad- 
justed empirically to provide the smoothest possible sky 
image. The quoted results are for a level of 5-a for sources 
detected at 24/im and 4-cr for undetected ones (see be- 
low), but they are not sensitive to modest adjustments in 
this level. In general, the resulting backgrounds have only 
weak structure and, where the sources are fairly bright 
in the infrared, the 24/im stacked image is similar to the 
point spread function of the instrument. We confirmed 
that stacking sources of known flux density gave consis- 
tent results. These behaviors validate the procedure. 

For this study, we used a total of 336 sources, selected 
similarly to those discussed above (but without screening 
for COMBO- 17 type classification). We divided the sam- 
ple into two groups. The first, hereafter the undetected 
group, includes 201 UV objects (60% of the total) for 
which visual inspection indicated no reliable 24/im detec- 
tion. The second (40%), hereafter the detected group, is 
the objects with evidence for an infrared detection; it was 
in turn divided into three equal subgroups according to 
NUV luminosity. As shown in Table 1, the two groups 
have very similar average NUV flux density, 1.89/iJy for 
the undetected and 1.63/iJy for the detected group. The 
average NUV luminosities are also similar, 1.3 x 10 44 erg/s 



and 1.1 x 10 44 ergs/s, respectively. However, the average 
infrared flux densities differ by a factor of ten: 13/iJy for 
the undetected group and 143/iJy for the detected one. 
There is no significant difference in average infrared flux 
density among the subgroups in the detected group. 

These results indicate that the LBGs divide into two 
classes. About 40% of them are infrared bright. The av- 
erage NUV flux density for this group is 1.63/iJy, so as 
measured in v¥ Ul the NUV and 24/im luminosities are 
similar. Since there is a substantial correction to total far 
infrared luminosity, the infrared component to the output 
from young stars is significant, probably accounting for the 
majority of the luminosity for these objects. The remain- 
ing 60% are infrared faint: i/F v is about ten times greater 
in the NUV than at 24/im, indicating that their outputs 
are dominated by the UV. The results of this paper apply 
to galaxies like those in the detected group only. 

To explore other possible differences between these 
classes, we computed the average B (i.e. rest-frame NUV) 
flux densities for the same two groups and three sub- 
groups. Although it is influenced by other factors, we take 
the ratio of B to NUV flux densities (or equivalently the 
NUV - B color) to be an indicator of the level of redden- 
ing, and the ratio of 24/im to B flux density to measure 
the relative portion of the luminosity from young stars 
emerging in the infrared compared with the UV. The re- 
sults are in Table 1. First, they demonstrate that all the 
galaxies in our selection fall in, or close, to the color-color 
LBG zone as adjusted from high z to z ~ 1 (see Giavalisco 
2002). There is a trend for LBGs with a high 24/B ratio 
to present a high B/NUV, which is consistent with the 
relation for the UV slope and the FIR/UV ratio found 
by Meurer et al. (2000) on a sample of local starburst 
galaxies. An analysis based on detections is required to 
check whether Meurer et al.'s law can be applied safely 
to those LBGs while we showed in the previous section 
that it provides dispersed SFRs for the detected sam- 
ple. Finally, the amount of dust attenuation for the un- 
detected group is very low (Apjjy ~ 0.5 — 0.6 for a mean 
LogLxoT ~ 10.5) which corresponds to LBGs with the 
lowest reddening found by Adelberger and Steidel (2000). 
This very low reddening is consistent with the very blue 
UV slope j3 » —2.4. If confirmed, this would mean that 
about half of the LBGs do not contain large amounts of 
dust. 

4. Conclusions 

We use multi-wavelength data in the CDFS to define the 
first large sample of Lyman Break Galaxies at z ~ 1; 
GALEX is used to observe the Lyman break. Redshifts 
are taken from spectra and from COMBO 17. Quantitative 
morphologies (Lauger et al. 2005a) are estimated from 
high spatial resolution images. Finally, dust attenuations 
and total luminosities are computed from SPITZER 
measurements at 24/im extrapolated to get the total IR 
luminosity. 

The main results of this analysis are: 
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Table 1. Stacking Analysis Results 





number 


Lf7V 


F(24^m) 


F(B) 


F(NUV) 


F(B)/F(NUV) 


F(24)/F(B) 






10 44 ergs/sec 


mjy 


mjy 


mjy 






All detected 


135 


1.11 


0.14 


0.0018 


0.0016 


1.11 


80 


low UV 


45 


0.67 


0.16 


0.0014 


0.0012 


1.21 


113 


middle UV 


45 


0.99 


0.16 


0.0017 


0.0014 


1.19 


95 


high UV 


45 


1.67 


0.14 


0.0023 


0.0024 


0.98 


59 


Undetected 


201 


1.33 


0.013 


0.0015 


0.0019 


0.79 


8.7 



1. We detect LBGs in the range 9.3 < LogLpuvlLG)] < 
11.0, i.e. well into the UVLG class as defined 
by Heckman et al. (2005). For the same objects, 
LogLm > HLq, which means that almost all the 
LBGs with a SPITZER confirmed detection are 
LIRGs (and 1 ULIRG) at z ~ 1. 

2. LBGs at z ~ 1 are mainly disk-dominated galaxies (75 
%) with a small contribution of interacting/merging 
galaxies (22 %) and a negligible (3 %) fraction are 
spheroids. The morphologies of our LBG sample are 
consistent with star-forming galaxies. 

3. The FUV dust attenuation appears to be anti- 
correlated with the observed FUV luminosity. Part of 
this correlation might be due to observational limits 
at 24 yon. However, non-detection of LBG with high 
Lpuv and high Apjjy cannot be explained by obser- 
vational biases. 

4. About 40% of our sample of LBGs is detected at 24/xm. 
These objects also show evidence for reddening in their 
UV continua. The remaining 60% of the sample are on 
average an order of magnitude less luminous in the in- 
frared compared with the rest frame near UV. Their 
UV continua also appear to be significantly less red- 
dened. 

5. Dust corrections of our IR-bright LBG sample com- 
puted via the j3 method are over-estimated by ~ 0.2 
mag. if we use Chary & Elbaz (2001) to compute 
Liu and by ~ 0.6 mag. if we use Takeuchi, Buat & 
Burgarella (2005). 

6. Dust attenuations estimated by the f3 method for such 
galaxies lead to overestimation of the SFRs at z ~ 1 
by a factor of 2 to 3, depending on the calibration of 
the 24/zm flux to Lju while IR-based SFRs are of the 
same order as SFRtot- 

7. By using the stacking method, we find that LBGs non 
detected at 24/xm seem to have very low dust attenu- 
ations. 
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